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Abstract. We present several variants of the Maupertuis principle, both on 
the exact and the non exact symplectic manifolds. 
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1. Introduction 

1.1. The principle of least action, or the principle of stationary action, says 
that the trajectories of a mechanical system can be obtained as extremals of a 
certain action functional. It is one of the basic tools in physics being applied both 
in classical and quantum setting. 

Consider a Lagrangian system (Q, L), where Q is a configuration space and 
L(q, q, t) is a Lagrangian, L : TQ xl->l. Let q = (qi, . . . , q n ) be local coordinates 
on Q. The motion of the system is described by the Euler-Lagrange equations 

ni , ddL dL 

(1.1) TiTT = ^—7 i = l,. ..,n. 

dt dqi dqi 

The solutions of the Euler-Lagrange equations are exactly the critical points 
of the action integral 



Sl(i) = / L(q,q, 

J a 



t)dt 



in a class of curves 7 : [a, b] — > Q with fixed endpoints 7(a) = qo, 7(6) = q\ (the 
Hamiltonian principle of least action (1834), e.g., see [28] ). 
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The Legendre transformation FL : TQ — > T*Q is defined by 




Pi = tt^, i = l,...,n, 



where £, r\ £ T q Q and (q±, . . . , q n ,Pi, ■ ■ • ,Pn) are canonical coordinates of the cotan- 
gent bundle T*Q. In order to have a Hamiltonian description of the dynamics (see 
the section below), we suppose that the Legendre transformation (|1.2[) is a diffco- 
morphism. The corresponding Lagrangian L is called hyperregular [21] . 

If the Lagrangian L does not depend on time then the equations (jl.ip possess 
the energy first integral 



In that case we have 

Theorem 1.1 (the Maupertuis principle). Suppose that h is a regular value of 
E. Among all curves q = 7(7") connecting two points qo and q\ and parametrized so 
that the energy has a fixed value E — h, the trajectory of the equations of dynamics 
(jl.ip is an extremal of the reduced action 



It is important to note that the interval [a, b], parametrizing the curve q — 7(7"), 
is not fixed and it can be different for different curves being compared, while the 
energy must be the same. 

Contrary to the Hamiltonian principle, the Maupertuis principle, or principle 
of stationary isoeneretic action determines the shape of a trajectory but not the 
time. In order to determine the time, we have to use the energy constant. 

Historically, a variant of Theorem 1.2 was the first variational approach to 
mechanics. It is attributed to Maupertuis (1744), Euler (1744) and Jacobi (1842), 
who gave an important geometric interpretation of the principle (see [28 ). 

1.2. The classical proofs of the Maupertuis principle can be found in 1281 1361 
[2]. In Serbian, see the second volume of Bilimovic's course in Theoretical mechanics 
[4], or Dragovic and Milinkovic's monograph |10j . 

Weinstein [34j and Novikov 25 formulated multi- valued variational principles 
that provided the study of the existence of periodic orbits on non exact symplectic 
manifolds. We feel a need to present these results, along with the classical ones, in 
a unified way. 

In the first part of the paper, we derive the principle of stationary isoenergetic 
action, both on the exact (Section 2) and the non-exact symplectic manifolds (Sec- 
tion 3). The variants of the Maupertuis principle presented in Section 3 are our 
small contribution to the subject. They slightly differ from the existing variational 
principles formulated either for closed trajectories, or formulated without imposing 
the constraint given by the energy. 

In the second part of the paper we point out a contact interpretation of the 
Maupertuis principle (Sections 4, 5). There, it is illustrated how some of the well 
known properties of the system of harmonic oscillators, the Kepler problem (Moser's 
regularization) and the Neumann system (relationship with a geodesic flow on an 



(1.3) 




(1.4) 
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ellipsoid), have natural descriptions within a framework of the contact geometry. 
We believe that one should expect other interesting relations between the contact 
structures and integrable systems as well. 

It is a great pleasure to dedicate this paper to Anton Bilimovic, since his work 
has fundamentally influenced the development of Serbian theoretical mechanics. 

2. Principle of stationary isoeneretic action in a phase space 

2.1. Hamiltonian equations. Let L(q,q,t) be a hyperregular Lagrangian. 
We can pass from velocities qi to the momenta pj by using the standard Legendre 
transformation p.2[) . In the coordinates (q,p) of the cotangent bundle T*Q, the 
equations of motion (jl.ip read: 

(21) dqi_ _dl£ dpi _ dH , 

dt dpi ' dt dqi ' 
where the Hamiltonian function H(q,p,t) is the Legendre transformation of L 

H(q,p,t) = E(q,q,t)\,j =¥L -i (q ^^ = ¥L(q,q,t) ■ q - L(q,q,t)\ i=WL -i^ Ptt y 

Let pdq = ^2iPidqi be the canonical 1-form and 

n 

lo = d(p dq) = dp A dq = dpi A dqi 

<=1 

the canonical symplectic form of the cotangent bundle T*Q. The system of equa- 
tions (|2.1[) is Hamiltonian, that is the vector field 

X H (q,p) = {dH/dpx, . . .,dH/dp n , -dH/dq u . . . , -dH/dq n ) 

can be defined by 

(2.2) i XH u{-)=w{X H , -) = -dH{-). 

2.2. Characteristic line bundles. More generally, consider a 2n-dimensional 
symplectic manifold P with a closed, non-degenerate 2-form lo. Let ff:Pxt^I 
be a smooth, in general time dependent, function. Consider the corresponding 
Hamiltonian equation 

(2.3) x = X H , 

where the Hamiltonain vector field Xn(x,t) is defined by (|2.2[) . 

If the Hamiltonian H does not depend on time, it is the first integral of the 
system. Let M be a regular connected component of the invariant variety H = h, 
which means dH\M 7^ 0. 

Since dH(£) = 0, £ e T X M, from (|2.2|) we see that Xh generates the symplectic 
orthogonal of T X M for all x £ M — the characteristic line bundle Cm of M. It is 
the kernel of the form u> restricted to M: 

C M = {£ G T X M I T X M) =0,i£ M}. 

Note that Cm is determined only by M and not by H . If F is another Hamilton- 
ian defining M, M C F _1 (c), d-F|M 7^ 0, then the restrictions of the Hamiltonian 
vector fields Xh and Xp to M are proportional. 

A variation of a curve 7 : [a, 6] — >• M is a mapping: L : [a, 6] x [0, e] — > M, 
such that 7(t) = L(t, 0), f G [a, 6]. Denote 7s (i) = L(i, s) and Sj(t) = ^| s=0 7 s (t) € 
T 7(t) M. 
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From Cartan's formula we get (e.g., see Griffits [12 ): 

Lemma 2.1. Let (M, a) be a manifold endowed with a 1-form a, 7 : [a, b] — > M 
be an immersed curve and T be a variation of 7 . The Lie derivative of the form 
T*a in the direction of d/ds at the points [a, b] x {0} is equal to 

Ld/dsF*ct\( t ,o) = l*{is 1 (t)da) + d^*(a(5j(t))). 

Theorem 2.1. Assume that the symplectic form uj is exact: uj — da. Let M 
be a regular component of the invariant hypersurface H~ l {h). The integral curves 
7 : [a, b] — !• M of the characteristic line bundle Cm are extremals of the (reduced) 
action functional 




in the class of variations j s (t) such that a (57(a)) = a(6~f(b)) = 0. 
The proof is a direct consequence of Lemma 2.1. We have 

(2.4) j- a) | s=0 = j a u } {6 1 {t)^{t))dt + a{8 1 {b))-a{8 1 {a)). 

The expression above is equal to zero for all variations 7 S (t) if and only if 7 is 
in the kernel of the form uj = da restricted to M . That is, j(t) is an integral curve 
of the line bundle Cm- 

2.3. Applying Theorem 2.1 to the symplectic space (T*Q, dp A dq) we obtain 
Poincare's formulation of the Maupertuis principle in a phase space [27 \. 

Theorem 2.2. // the Hamiltonian function H = H(q,p) does not depend on 
time, then the phase trajectories of the canonical equations (|2.1j) lying on the regular 
connected component M of the surface {H{q,p) = h} are extremals of the reduced 
action 

(2.5) A(r() = f pdq 

in the class of curves 7 lying on M and connecting the subspaces T* g Q and T* t Q. 

Note that Theorem 1.1 follows from Theorem 2.2 (e.g., see Arnold [2]). Suppose 
that the Hamiltonian system (|2.1j) is a Legendre transformation of the Lagrangian 
system (|1.1|) . The main observation is that if 7(7-) is a configuration space curve 
parametrized such that £^7 ,dj j 'dr) = h, then the lifted curve 7 = ¥L(~f, d'y/dr) 
lyes on M and the reduced actions (|1.4[) and f|2.5[) for 7 and 7 are equal: ^(7) = 
^(7) (see Fig. 1). 

2.4. Jacobi's metric. Consider a natural mechanical system on Q defined 
by the Lagrangian function: 

(2.6) L(q, q )=T + B-V=^Y, K M + E B <* " 

ij i 

Here ds 2 = ■ Kijdqidqj is a Riemannian metric on Q, V(q) is a potential function 
and 6 = £\ Bidqi is a 1-form defining a gyroscopic (or magnetic) field a = d9 (see 
Section 3). 
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Figure 1. 



The energy of the system (|1.3p is the sum of the kinetic and the potential 
energy 



In the region of the configuration space Qh where V(q) < h, we can define the 
Jacobi metric 



The following version of the Maupertuis principle for Lagrangians of the form 
(|2.6p is well known (e.g., see Kozlov [19 ). 

Theorem 2.3. Among all curves q = 7(t) connecting the points qo,q\ G Qh 
and parametrized so that the energy has a fixed value E = h, the trajectory of the 
equations of dynamics (|1.1[) with Lagrangian (12.61) is an extremal of the integral 



In particular, if there are no gyroscopic forces, the trajectories of the system within 
Qh, up to reparametrization, are geodesic lines of the Jacobi metric ds 2 j. 

Indeed, in order to guarantee a fixed value of the energy 



E(q, q) = T+V=-J2 KiMi + V(q). 



(2.7) 



dsj = 2(h - V(q)) ds 2 = 2(h - V(q)) ^ K^dqdq^ 



(2.8) 




7 



ij 
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the parameter r of the curve q — j(t) must be proportional to the length dr — 
ds/y/2(h-V). Therefore 




Remark 2.1. The variational principle stated in Theorem 12.31 is used in the 
study of periodic trajectories of natural mechanical systems with exact magnetic 
fields (see [31] and references therein) . Note also that the Maupertuis principle for 
a configuration space Q being a Banach space can be found in |21L I30j . 

2.5. The Hamiltonian principle of least action. Consider a Poincare- 
Cartan 1-form pdq — Hdt on the extended phase space T*Q x WL(q,p, t), where 
H : T*Q x M. — >• K is a Hamiltonian function. The phase trajectories of the 
canonical equations (I2.1j) are extremals of the action 

(2.9) A H (rf) = f pdq -Hdt 

in the class of curves j(t) — (q(t) , p(t) , t) connecting the subspaces T* Q Q x {to} and 
Tg t Q x {ti} (Poincare's modification of the Hamiltonian principle of least action 
[27)). Namely, a vector (£,1), £ € T {q ^ p) (T*Q) belongs to kevd(pdq - Hdt) at 
(q,p,t) if and only if £ = X H (q,p, t) (see [1 [21]). 

Obviously, we can replace (T*Q,dp A dq) by an arbitrary exact symplectic 
manifold (P, co = da). In particular, if we consider the action ^#(7) = a — Hdt 
on the free loop space fi(P) = ^(S 1 ,?), S 1 = R/Z of P and P is 1-periodic in 
t- variable, then the critical points of Ah are 1-periodic orbits of the equation (|2.3p . 

For a given time-independent Hamiltonian H : P — ¥ K with a regular level 
set P _1 (/i), the periodic orbits having all positive periods and energy h can be 
obtained by the use of modified action: 

(2.10) A Bth (rf,X)= [ a(j)dt-\ [ (H( 7 (t))-h)dt, 

Jo Jo 

defined on the space fi(P) x R + (see [29l [34]). The critical points (7, A) of A H ,h 
correspond to A-periodic orbits x(t) = 7(t/A) that lie on the energy hypersurface 
H^ 1 (h). Moreover, Weinstein defined actions Ah and Ah,k when the symplectic 
form is not exact as well [35] . 

The Lagrangian analogue of the functional (12 . 10[) is 

S L ,h(l> A ) = / XL h-, i/^)dt + Xh, 7 € n(Q), A > 
Jo 

(see [9]). The pair (7, A) is a critical point of Sl,h if and only if q(t) = j(t/X) is a 
A-periodic solution of the Euler-Lagrange equation (|1.1[) with energy h. 

Variational principles related to the action (|2.9p . which arise by a reduction 
process are given in [8]. 
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3. The Maupertuis principle on non-exact symplectic manifolds 

3.1. Magnetic flows. Consider a natural mechanical system given by the 
Lagrangian function (|2.6p . After the Legendre transformation, it takes the form 
(|2.ip with the Hamiltonian function 

(3.1) H(q,p) = l(p-6,p-6)+ V(q) = ^ " - *i) + V(q), 

ij 

where K l i is the inverse of the metric tensor Kij . 
The transformation 

(3.2) Tg:(q,p)^(q,p-0) 

is a symplectomorphism between (T*Q, dpAdq) and a "twisted" cotangent bundle 
(T*Q, dp A dq + 7r*cr), where 7r : T*Q — > Q is the natural projection and a = dQ. 

In new coordinates, also denoted by (q,p), the Hamiltonian (|3.1[) takes the 
usual form, the sum of the kinetic and the potential energy: 

H(q,p) = l(p,p) + V(q) = lj2 Kij P^ +V (ti> 

ij 

while the equations of motion take "non-canonical" form: 

dq* _dH d Pi __dH ^ OH 
[ > dt ~ d Pi ' dt ~ ihr ' ^ '<)Pj' 

where a = Y^,i<i<j^ n Fij(l)dqi A dqj. The equations are Hamiltonian with respect 
to the symplectic form oj = dp A dq + n*a. 

One can consider the system (|3.3p associated to a non-exact 2-form a as well (for 
example, the motion of a particle in a magnetic monopole field |21j ). In this case, 
the Lagrangian (|2.6I) is defined only locally. Nevertheless, it is very interesting that 
the Hamiltonian ( Weinstein |34] and Tuynman |33j ) and the Maupertuis principles 
(Novikov [26 ) of least action can be still defined. 

3.2. Multivalued reduced action. Let (P, oj) be a non exact symplectic 
manifold and let M = _ff _1 (/i) be a regular isoenergetic hypersurface. The main 
observation concerning the Maupertuis principle can be stated as follows (see |26L 
H9] for the reduced action l|2.8pV 

Let U C P be a region where uj is exact and let uj — da\ — dct2- Consider a 
variation J s (t) = T(t, s), t £ [0, 1], s € [0, e] with fixed endpoints of a curve 7 lying 
in M n U. Then 



0(i— cti = — T*dai = — I T*da2 = I aa — / oli- 

'To ^[0,l]x[0,e] J[0,l]x[0,c] J~< € J~<o 



Therefore 



oli — I oli = / a 2 — a-2 

7o ■'To 



and, although J ai depends on the form a^, the derivative 



YJ s =°J q a (%) dt = J ^l{t)Mt))dt 

does not depend on c^, i = 1,2. One can define an appropriate multi- valued 
functional on a space of paths with fixed endpoints, such that an extremal (if exist) 
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is exactly the integral curve of the characteristic foliation on M. However, as in the 
case of the symplectic homology (see \13\). the situation simplifies in the aspherical 
case which is considered below. 

3.3. Aspherical symplectic manifolds. The symplecic manifold {P,oj) is 
aspherical if uj vanishes on %2(P)- Ofcourse, if uj is exact or ^(P) — 0, then (P, ui) 
is aspherical. 

Consider the equation (j2.3[) , where H does not depend on time. Let M be 
a regular component of -ff -1 (/i) and c : [0, 1] — > P be an immersed curve with 
endpoints xq = c(0) € M and x\ = c(l) € M. Define f2^(a;o,xi) as the space of 
regular paths that are homotopic to 7 in P: 

Q h c (x , n) = {7 : [0, 1] -> M I 7 (0) = x , 7(1) = x 1; 7 (i) ^ 0, t € [0, 1], 7 ~ P c}. 

The space of all regular paths connecting xq and x\ and laying in M is the 
union 

n h (x ,x 1 ) = \Jn%(x ,x 1 ), 

c 

where we take representatives c for all non-homotopic paths (in P) connecting xq 
and X\. 

If we suppose that (P, uj) is simplectically aspherical then we can define a 
single-valued reduced action: 

(3.4) A : Q h (x ,x 1 ) -► E, A( 7 )|^ (x0jKl) = / />, 

where 

D = {z|z e C, |z| < 1} 
is the unit disk, / 7 : D — > P is an arbitrary mapping that is smooth for \z\ < 1, 
continuous on D and 7 (t) = / 7 (exp( v / — l7r£)), c(t) — / 7 (exp(V~ l7r(2— <)), t € [0,1]. 
That is, /(-D) is a surface with the boundary dD — 7 • c _1 . 




Figure 2. 

Since 7~pc we can always find a mapping / with required properties. From 
UJ \n2(P ) = 0; the value ^(7) does not depend on the choice of /. 

Theorem 3.1. The integral curves 7 : [0,1] — > M of the characteristic line 
bundle Cm that connect xq and x\ are extremals of the reduced action (13.41) . 
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Proof. Consider a variation %(t) — T(t,s), t £ [0, 1], s £ [0, e] of 7 lying in 
M. By using u\„ a ( P ) = 0, we get 



A( 7e ) - A( 7 ) = / (/> - /» = - 
Thus, as above, 



[0,l]x[0,s] 



d 



, „=o^(7 s ) - / w(*7(t),7(t))dt, 

is zero for all variations 7s (t) if and only if the velocity vector field 7 (i) is a section 
of ker u)\m- O 

3.4. A torus valued reduced action. Tuynaman proposed a torus-valued 
action, such that the multi-valued Poincare action (12.91) can be seen as a compo- 
sition of a multi- valued function on a torus and a torus- valued action |33j . In 
this subsection we follow Tuynman's construction |33j in order to formulate the 
principle of stationary isoenergetic action. 

Consider a manifold P with a symplectic 2-form 

n 
a=l 

where f3 a are 2-forms, representing integrals cohomology classes. We take decom- 
position with minimal n. Then the parameters fi a are independent over Q, in 
particular jj, = fj,% + • • • + /i„ 7^ 0. To lo we associate the 1-form 



A = ^dy a 



a=l 



on a torus T" = {(cxp(-\/^Iy 1 ), . . . , exp(-\/^Ty™)}. It can be consider as a differ- 
ential of a multi- valued function A on T": A = dA. Also, for a = 1, . . . , n, let us 
define principal ^-bundles 

S 1 — ► ^ 

(3.5) 



having the connections 6 a with the curvature forms j3 a (see Kobayashi (18] ). 

Let 7 (t), t £ [to,t\\ be a piece-wise smooth, closed curve on P. Recall, a 
piece- wise smooth curve j a (t) C Y a is a horizontal lift of 7 if ,o a o 7 a (t) = 7(i) and 
® a (lil a (t)) — 0' whenever the velocity vector is defined. The holonomy Hol a ( 7 ) is 
an element g £ S 1 , such that g ■ 7 a (io) = l a (ti)- 

Lemma 3.1 ([33]). Let 7s (i) = F(t,s) be a variation of 7 : [0,1] — > P with 
fixed endpoints and let c : [0, 1] — > P be an arbitrary curve connecting xq = 7(0) 
and X\ = 7(1). We have a family of closed orbits % = j s ■ c . The derivative of 
Hol a (7 s ) is given by: 



ds ./,, ''//'•" 
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Consider the equation (|2.3[) . where dH/dt = 0. Let M be a regular component 
of H~ x (h) and Q h (xq, x\) be a space of regular paths 7 : [0, 1] — ¥ M that connect 
points X\ and £2- 

For every 7 G ^(xo, xi) define a pease- wise smooth, closed path 7 = 7 • cT 1 : 
[0,2] ->■ M, where c G Q h (x ,ai) is fixed. We call 

A Tn : £l h {xo,Xx) — ¥ T™, 7^ (Hoi 1 (7),..., Hoi" (7)) 

a toras valued reduced action. From Lemma I3JJ we have: 



AQUr»( 7 .)Uo) = X>«j(V(i 



(M/UM/)W/ = / «(7(i),« 7 (t))dt. 



Whence, we obtain the following principle of stationary isoenergetic action on 
the non-exact symplectic manifolds. 

Theorem 3.2. ^4 curve 7 G f2' l (xo, £1) is an integral curve of the characteristic 
line bundle Cm if and only if 

^ (Ao^y | s=0 = A f Aa t »( 7s )| s=0 ) = 

/or a/Z variations 7 S G Q (xo, Xi). 

For the completeness of the exposition we include: 

Proof of Lemma 13.11 In local trivializations 

/5 _1 (^() = UiX S 1 (x i ,y i mod 2%), 

we have local connection 1-forms cti on Ui such that 6 = on + dyi (the index a 
is omitted). The transition functions between fiber coordinates and connection 
1-forms are given by 

(3.6) y =y t + gi 3 {x), 04 = aj + dg^, g^ : U l n Uj -¥ S 1 . 

On the other hand, the curvature 2- form is invariant: /3 = dcti = dctj. 
Suppose 7a([io,ii]) C Ui, s G [0,e]. The local expression for % reads 

7-(*) = (7-(t).»i(*.*)). e *(^(*)) =0 a .(7(*))+^(M) = 0. 

Therefore 

(3.7) Vi(h,s) = Ui(to,s) - I ai(j s (t))dt. 

J to 

By taking the differential of (|3.7[) at s = and applying (|2.4[> we get 

(3.8) %(ti) + ai(*y(ti)) = S yi (t ) + a^fo)) + f * P(i(t), ^{t))dt, 

Jtg 

where o";y(t) = s)U=o, fry(t) = ^7«(*)U=o- 

Now, assume t < t < *i < 7a ([*0)*i]) C £/* and 7 s ([t 0) *i]) C t/j. The 
transformations (|3.6I) imply 

(3.9) *tfi(*) + o<(*7(*)) = <%(*) +a J -(<J 7 (t)) ) * G [4,t x ]. 
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By combining p.8[) and p.9[) , it follows 
(3.10) 5 yj (t\) + ai (*7(ti)) = «w(t ) + Oi(<y7(*o)) + / 1 /3(7(*)^7(*))^- 



Let 7s = 7 S • c" 1 : [0, 2] — > M and let U\, . . . , Ui be local charts, such that 
%([U-i,ti]) C U h = t Q < t x < ■■■ <t k = Kt k+ i < ■■■ < t t = 2, s€[0,e]. 
From the relation (jXTUf and (57(0) = £7(1) = = Sj(t) = 0, t € [1, 2], we get 

6y k (l)-6yx(0)= [ /?( 7 (i), M*))*, %(2) - <%(1) = 0. 
Jo 

We can suppose that the horizontal lifts of all curves start from the same point in 
Y. Then 6yx{0) — 0. This proves the statement. □ 

3.5. Reduced action for magnetic flows. Let us return to the magnetic 
equations p.3p . where H(q,p) is an arbitrary smooth function and a is not exact. 
Let M be a regular component of H(q : p)~ 1 (h) and let ir : T*Q — > Q be the natural 
projection. 

As in Theorem 2.2, we not need to fix endpoints in the fiber directions. Consider 
a class of regular curves 7 lying on M and connecting the subspaces T*Q and T* ± Q, 
such that the projection 77(7) is homotopic to c: 

(90, gi) = {7 : [0, 1] ->• M I tt( 7 (0)) = go, tt( 7 (1)) = nAl) ~ 4, 
and a class of all regular paths connecting T* Q and T* ± Q and laying in M: 

n h (q ,qx)=\Jn h c (q ,qx), 

c 

where we take representatives c : [0, 1] — > Q for all non- homotopic paths connecting 
g and 91 . 

Theorem 3.3. Assume (t\ V2 (q) = 0. The phase trajectories of the magnetic 
equations (|3.3|) m i/ie c?ass of curves O,^(qo,qx) are extremals of the reduced action 

A:Q h (q ,qx)^R, A( 7 )\ nH<]0tqi) = [ P dq+ [ f*a, 

J 7 J D 

where / 7 : D — » Q is smooth for \z\ < 1, continuous on D and 

*(7(t)) = / 7 (exp(V Z T 7 rt)), c(t) - / 7 (exp(^l7r(2 - <)), t e [0, 1]. 

If t| 7I - 2 (q) ^ 0, we can use a combination of the usual reduced action and a 
torus valued action with respect to the form a. Suppose 

n 
o=l 

where (3 a are 2-forms, representing integrals cohomology classes in Q. We take 
decomposition with minimal n. As above, to a we associate principal S' 1 -bundles 
L a over Q having the connections 6 a with curvature forms /3 a , a = 1, . . . , n. 

Let us fix c : [0, 1] —> Q, c(0) = go, c(l) = gi. For every 7 e fi' l (go,gi), we 
associate a pease- wise smooth, closed path 7 = 71(7) • c _1 : [0, 2] — > Q. Define 

B rn : £l h (qo, qx) — > T n , 7 .— ► (Hoi 1 (7), • ■ • . Hol n ( 7 )), 
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where Hol a is the holonomy of the bundle L a — > Q. Let v — X)"=i l JL ady a be a 
1-form on T™, considered as a differential of a multi- valued function T: dT = v. 

Theorem 3.4. A curve 7 £ fl h (qo,qi) is an integral curve of the characteristic 
line bundle Cm if o,nd only if 



for all variations 7 S G £l h (qo, qi). 

Remark 3.1. For various approaches to the existence problem of closed mag- 
netic orbits, see [9], 1311 and references therein. Integrable magnetic geodesic flows 
on homogeneous spaces can be found in [6]. 



4.1. A contact form a on a (2n+l)-dimcnsional manifold M is a Pfaffian form 
satisfying a A (da) n ^ 0. By a contact manifold (M, %) we mean a connected (2n + 
l)-dimensional manifold M equipped with a nonintcgrable contact (or horizontal) 
distribution locally defined by a contact form: H\u = kera\u> U is an open set 
in M [20 . A contact manifold (M,1-C) is co-oriented (or stricly) contact if H is 
defined by a global contact form a. For a given contact form a, the Reeb vector 
field Z is a vector field uniquely defined by 



4.2. In studying the existence problem of closed Hamiltonian trajectories on 
a fixed isoenergetic surface, Weinstein introduced the following concept [35] • An 
orientable hypersurface M of a symplectic manifold (P, w) is of contact type if there 
exist a 1-form a on M satisfying 



where j : M P is the inclusion. If (M, a) is of contact type, since C — ker ljm, 
the kernel of a 



is a (2n — 2)-dimensional nonintegrable distribution on which da = lo is nondegen- 
erate. Consequently, a A da 71-1 is a volume form on M and (M, "H) is a co-oriented 
contact manifold. 

Now, let (P, cj = da) be an exact symplectic manifold. Consider a regular 
component M of an isoenergetic surface (i? does not depend on time). If 

o.{Xh)\m 7^ then M is of contact type. We say that M is of contact type with 
respect to a. 

If M is of contact type with respect to a then a has no zeros in some open 
neighborhood of M. Contrary, suppose that an 1-form a has no zeros in some open 
neighborhood of M. Then, from the nondegeneracy of lj, there exists a unique 
vector field E such that 

(4.1) isuj = a. 

The vector field E has no zeros. From Cartan's formula, the condition i^cj = a is 
equivalent to L^w = uj, i.e., E is the Liouville vector field of lo. We have (e.g., see 
Libermann and Marie [20] ) : 




4. Isoenergetic hypersurfaces of contact type 



i z a = 1, 



izda = 0. 



n = {£ G T X M I a(f) =0,ie M} 
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Lemma 4.1. A regular connected component M of an isoenergetic surface -ff -1 (/i) 
is of contact type with respect to a if and only if the Liouville vector field defined 
by (|4.1|) is transverse to M . 



Proof. Since iE^ n — naAda n 1 , the kernel of a A da™ 1 is the vector bundle 
generated by E. Therefore a A da n \m is a volume form on M at x if and only if 
E{x) i T X M. □ 

Let M be of contact type with respect to a and let Z be the corresponding 
Reeb vector field on M: 

izda\ M = 0, a(Z) = 1. 

Since Z is a section of kerrfa|A/, it is proportional to Xh\m- Z = MXh\mi 
Af 7^ 0. Consequently, the flow of Z can be seen as a flow of Xh\m after a time 
reparametrization dt = Afdr: 

(4.2) ^ = dxdt = ^ ^ x € M. 
dr dt dr 

Alternatively, we can change the Hamiltonian H. Extend M to a neighborhood of 
M. Then 

(4.3) X mH _ h) {x)=M(x)X H (x), xeM. 

Based on the observations (|4.2[) . (|4.3p . we have the following statement. 
Lemma 4.2. The function 

H-h 
H °-E{H) 

has M as an invariant surface and the Hamiltonian vector field Xh \m is equal to 
the Reeb field Z . If p is any smooth function of a real variable, such that p'(X) = 1, 
then p(Hq+X) has the same property. In particular, for p{x) = —l/(4x), A = —1/2, 
we get 

(4-4) Hj = Ah-AH H +2E{Hy Hj]m = 5' Z = X *>\™' 

PROOF. According to (|2.2|l . (|4.ip . we have 

(4.5) a{X F ) = uj(E, X F ) = dF(E) = E(F), F G C°°(P). 

Thus, Z = X h /E{H)\m, i.e., Af = l/E(H). It is clear that H \ M = 0, while 

implies Z = A# | Af . 
Let p is a smooth function, such that p'(X) = 1. Then /o(Ho + A)|m = p(A) and 
£(p(ff + A))|M=p'(A)£(!/ ) = l. ' □ 

4.3. Exact magnetic flows. Consider a natural mechanical system given by 
the Hamiltonian function (|3.1I) . The canonical 1-form pdq is different from zero 
outside the zero section {p = 0}, where we have the standard Liouville vector field 
E = T,iPid/d Pi on T*Q. 

Since E(H) = (p,p — 6), a regular hypersurface Mh = H^ 1 (h) is of contact 
type with respect to pdq within a region 

M , h - {{p - 9,p - 9) + 2V(q) = 2h, (p,p- 9)^0} 

= {{ P -9,p-9)+2V(q) = 2h, (p,p) j^2V + (6, 6) - 2h} C T*Q h . 
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Note that the equation (p,p — 9) = 0\ q , 9 q ^ 0, defines an ellipsoid in T*Q. 
Assume 

K = max ( V(q) + -{9, 9) ) < oo 

q£Q \ 2 J 

(for example, h* exists if Q is compact). Then regular hypersurfaces M — H (h), 
for h > ft,*, are of contact type with respect to pdq. The function (|4.4p has the form 

Hj(q,p) = (p- e,p)/(i(h - V{q)) + 2(9,p)). 

In particular, if 9 = 0, Hj is the Hamiltonian function of the geodesic flow of 
Jacobi's metric (|2.7[) and Mh is the corresponding co-sphere bundle over Q. 

Remark 4.1. The function Af in the time reparametrization (|4.2I) equals Af = 
l/E(H), E(H) = (p,p - 0) = 2(h - F(<z))|m. That is, 

dt = dr/2{h - V), 

which agrees with Corollary 2.2 (where the time parameter dt of the original system 
is denoted by dr, dr = ds/y / 2(h — V) = dsj/2(h — V), and dsj is the natural 
parameter of Jacobi'metric). 

5. Examples: contact flows and integrable systems 

5.1. Harmonic oscillators. Consider the simplest integrable system - the 
system of n independent harmonic oscillators defined by the Hamiltonian function 

(5.1) H = J2 F *> Fi = \{ad + hp1)i l,...,n, 

i 

in the standard symplectic linear space M. 2n (q,p). Here we suppose that the prod- 
ucts ai&i, i = 1, . . . , n are positive. 

By the use of the first integrals Fi — Ci, a generic solution of the equations 

(5.2) fa = biPi, pi = -cnqi, i = l,...,n 

can be written in the form qi{t) = cos + </j°), Pi(t) — — y^^sin {ojit + ipf) , 
0Ji — y/dibi, where (p® £ [0, 2ir) are determined from the initial conditions. Assume 
A k = a ri+ ... +Tk _ 1+1 = ■ ■ ■ = a,. lH — hrfc , 

Bk = &nH hr-fc-i+1 = • • • = \- rh , 

1 < k < s, ri H h r s = n, ro = 

and that the frequencies \fA\B[, \J A2B2, . . . , ^/ A S B S are independent over Q. 

Due to the U(r\) x ••• x {J(r s )-symmetry, the system (|5.2I) has additional 
Noether integrals 

Fij = - \ 1.(1,(1, + I'll', l>,. = '/./'. - /','/-• 

J"H h r fc _i + 1 < i < j < n H h r fc , fc = 1, . . . , s, 

implying the non-commutative integrability of the system |25l 122) . Generic trajec- 
tories fill up densely invariant s-dimensional invariant isotropic tori generated by 
the Hamiltonian vector fields of integrals 

H\= F\ + ■ ■ ■ + F ri , . . . , H s = F ri -\ hr s _i+l + • • • + -FriH H's • 

The quadric Mh = i/ _1 (ft), h ^ is of contact type with respect to the 
canonical 1-form pdq outside p = 0, where we have a well defined Jacobi's metric. 
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However, if instead of pdq, we take 

n 1 " 1 

(5-3) a = ^2p l dq i - ^(^PiQi) = , $^ p<d * - Qidpt, 

i—1 i—1 i 

then da — d{pdq) — dp A dq and the only zero of a is at the origin 0. The 
corresponding Liouville vector field is 

Since E(H) = h\M h , the quadric Mh is of contact type with respect to a and 
the Reeb flow on Mh is Z = h~ 1 Xfi\M h - 

The above construction provides natural examples of contact structures on 
quadrics within R 2 ™ having the integrable Reeb flows with s-dimensional invariant 
tori, for any s — l,...,n. The case s — n corresponds to contact commutative 
integrability introduced by Banyaga and Molino [3] (see also |16L [7]), while for 
s < n we have contact noncommutative integrability recently proposed in [14] . 

By taking all parameters to be positive (oj, bi > 0, i = 1, . . . , n), after rescaling 
of Mh to a sphere S* 2 ™ -1 , we get i^-contact structures on a sphere iS 2 ™ -1 given by 
Yamazaki (see Example 2.3 in [37] ). In particular, for a% = ai = • • • = a n = b n = 1 
we have the standard contact structure on a sphere 5 n_1 = H~ 1 (l/2) with the 
Reeb flow which defines the Hopf fibration (e.g., see [20] ). 

Remark 5.1. A modification of the canonical form pdq given by (|5.3p can be 
applied for starshaped hypersurfaces in R 2 ™. More generally, consider a regular 
isoenergetic hypersurface Mh = H~ x in (T*Q(q,p),dp A dq). It is of contact type 
if there exist a closed 1-form <p on Mh such that pdq{Xn\M h ) + ^p{XH\M h ) ^ 
0. Assume Mh is compact. Then the required 1-form ip exists if and only if 
J M pdq(Xn)dii 7^ for every invariant probability measure \i with zero homology 
(see Apendix B in [§]). In particular, for a compact regular energy surface Mh = 
ff _1 (/i) in the standard symplectic linear space (R 2n (q,p),dq A dq) we have the 
following sufficient conditions. Suppose: 

(i) pdq(X H )>0, forp^O, (q,p)£M; 

(ii) if M n {p = 0} ^ 0, then ^H(q, 0) ^ at the points (q, 0) G M. 
Then Mh is of contact type with respect to 

n ( n d \ 

a = '^p l dq l - ed I ^ft— H(q,0) J , 

i—1 \i— 1 / 

for a certain parameter e (see [13] ). 

5.2. The regularization of Kepler's problem. The motion of a particle 
in the central potential filed is described by the Hamiltonian function 

H:Rl n =R 2n ^{q = 0}^R, H(q,p) = M- - SL 

where (•, •) is the Euclidean scalar product in R™. Moser's regularization of Kepler's 
problem (see [23] ) can be interpreted in contact terms as follows. 

Let Mh = {H = h} C R 2rl be an isoenergetic hypersurface. Let us interchange 
the roll of q and p and consider the form a — — X)™=i Qidpi and the associated 
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Liouville vector field 

^ d 



1=1 

Since E(H) — "i/\q\, Mh is of contact type with respect to a. According to 
Lemma |4~21 the Reeb flow on Mh can be seen as a Hamiltonian flow of 

Hv = {\p\ 2 -2h)\q\/2 1 -l. 

In order to get a smooth Hamiltonian we can take F = (Hq + 1) 2 /2 (see Lemma 
E2): 



2 



Then F\M h = \, Z = Xp\M h and, moreover, Xp is defined on the whole R 2 ™. 
Assume h < 0. The Hamiltonian F(q,p) can be interpreted as a geodesic flow 
of the metric proportional to 

2 _ d Pj + • • • + d Pl 
h (2/1 -bl 2 ) 2 ' 

It represents the round sphere metric obtained by a stereographic projection (see 
Moser [23] ). Thus, for h < 0, there exist a compact contact manifold Mh = 
Mh U 5™ (a co-sphere bundle over S n ) with a Reeb vector field Z, which is a 
smooth extension of Z. In particular, for n = 2, M/j = RP 3 . On RP 3 we have a 
standard contact structure, obtained from the standard contact structure on S 3 via 
antipodal mapping. 



Note that for h > 0, the metric ds^ is defined within the ball of radius y/2h 
and represents Poincare's model of the Lobachevsky space. 

The contact regularization of the restricted 3-body problem is given in pQ . 

5.3. The Maupertuis principle and geodesic flows on a sphere. It is 

well known that the standard metric on a rotational surface and on an ellipsoid 
have the geodesic flows integrable by means of an integral polynomial in momenta 
of the first (Clairaut) and the second degree (Jacobi) [2j. A natural question is 
the existence of a metric on a sphere S 2 with polynomial integral which can not 
be reduced to linear or quadratic one. The first examples are given in [5] . Namely, 
the motion of a rigid body about a fixed point in the presence of the gravitation 
field admits <SO(2)-reduction (rotations about the direction of gravitational field) 
to a natural mechanical system on S 2 . Starting from the integrable Kovalevskaya 
and Goryachev-Chaplygin cases and taking the corresponding Jacobi's metrics, we 
get the metrics with additional integrals of 4-th and 3-th degrees, respectively. 

We proceed with a celebrated Neumann system. The Neumann system de- 
scribes the motion of a particle on a sphere (q, q) = 1 with respect to the quadratic 
potential V(q) — h(Aq,q), A = diag(oi, . . . , a n ) (we assume that A is positive 
definite). The Hamiltonian of the system is: 

(5.4) H N (q,p) = -(p,p) + -{Ax,x). 

Here, the cotangent bundle of a sphere j^g"- 1 i s realized as a submanifold P 
of R 2 ™ given by the constraints 

(5.5) F 1 = (q,q) = l, F 2 = (q,p)=0. 
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The canonical symplectic form on P = T*S n ~ 1 is a restriction of the standard 
symplectic form dp A dq to P. Let H : M. 2n — > R. The Hamiltonian vector field 
Xh\p reads 

^h(<7,p)|p = x H(q,p) - XiX Fl (q,p) - X 2 X F2 (q,p), (q,p) e P, 

where the Lagrange multipliers are determined from the condition that Xjj\p is 
tangent to P (e.g., see [24]). 

There is a well known Knorrer's correspondence between the trajectories q(t) 
of the Neumann system (|5.4|) restricted to the zero level set of the integral 

(5.6) H(q,p) = ~ ((A- 1 q,q)(A- 1 p,p) - (A^q.^p) 2 - (A^q)) . 

and the geodesic lines on an ellipsoid E™ -1 = {x S W 1 \ (x, Ax) = 1} by the use of 
a time reparametrization and the Gauss mapping q — Ax/\Ax\ |17j . 

Recently, by using optimal control techniques, Jurdjevic obtain a similar state- 
ment for the flow of the system defined by the Hamiltonian (|5.6[) 15 . 

We give the interpretation of Jurdjevic's time change by the use of Maupertuis 
principle. Since the potential V(q) = —^(A~ 1 q,q) is negative, the isoenergetic 
surface 

(5.7) M = {H\ P = 0}cP = T*S n ~ 1 

is of contact type with respect to pdq\p. The Reeb vector field Z equals to the 
Hamilonian vector field of 

(5-8) Hj = — 1 {{A-\,q){A- l p,p) (A-^pf) \ P 

4(A 1 q,q) 

(the Hamiltonian of the corresponding Jacobi's metric). 

The Legendre transformation of a function of the form (|5.8p in the presence of 
constraints (I5.5[) is given in (see Theorem 2 and interchange the role of the 
tangent and cotangent bundles of a sphere). As a result, we obtain the Lagrangian 
function 

L(q,q) = 2<^g,9)|s»-i • 

Remarkably, after the linear coordinate transformation x = \f~Aq, L(q,q) be- 
comes the Lagrangian L(x,x) = ^(x,x) of the standard metric on the ellipsoid 

K' 1 = {xe M™ | (A~ 1 x,x) = 1}. 

Recall that the Reeb flow on Mo can be seen as a time reparametrization of the 
original Hamiltonian flow (see Remark 4.f). We can summarize the consideration 
above in the following statement. 

Proposition 5.1 ([H]). Under the time substitution 

dt = dT/2(A~ 1 q,q) 

and the linear transformation x — \f~Aq, the q-components of the trajectories of the 
system defined by the Hamiltonian function (I5.6[) that lye on the zero energy level 
(|5.7p . become geodesic lines of the standard metric on the ellipsoid E^~ x . 

Further interesting examples of transformations related to the Maupertuis Prin- 
ciple, which map a given integrable system into another one are given in [32 . 
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